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ABSTRACT 

Available X-ray data are collected and organized concerning the iron and gas 
content of galaxy clusters and groups, together with the optical luminosity, mass 
and iron abundance of cluster galaxies. Moving from such a restricted number of 
cluster parameters several astrophysical inferences are drawn. These include the 
evidence for rich clusters having evolved without much baryon exchange with 
their surrondings, and having experienced very similar star formation histories. 
Groups are much gas-poor compared to clusters, and appear instead to have 
shed a major fraction of their original cosmic share of baryons, which indicates 
that galaxy clusters cannot have formed by assembling groups similar to the 
present day ones. It is argued that this favors low-G universes, in which the 
growth of rich clusters is virtually complete at high redshifts. It is also argued 
that elemental abundances in clusters are nearly solar, which is consistent with 
a similar proportion of supernovae of Type la and Type II having enriched 
both the solar neghborhood as well clusters as a whole. Much of the iron in 
clusters appears to reside in the intracluster medium rather than inside galaxies, 
the precise ratio being a function of the Hubble constant. It appears that the 
baryon to star conversion in clusters has been nearly as efficient as currently 
adopted for the universe as a whole. Yet the metallicity of the clusters is ~ 5 
times higher than the global metallicity adopted for the nearby universe. It is 
concluded that the intergalactic medium should have a metallicity ~ 1/3 solar 
if stellar nucleosynthesis has proceeded in stars within held galaxies with the 
same efficiency as in stars within clusters of galaxies. 

Subject headings: galaxies: abundances - galaxies: clusters: general - galaxies: 
formation - galaxies: intergalactic medium 



1. Introduction 


The X-ray observations of clusters of galaxies have revealed the presence of large 
amounts of iron and other heavy elements in the intracluster medium (ICM) (Mitchell, Ives, 
& Culhane 1975; Serlemitsos et ah 1976; Mushotzky et ah 1996), thus providing direct 
evidence that gas contaminated by nucleosynthesis processes has been lost by galaxies in 
the course of their evolution. The access to the ICM elemental abundances has offered the 
opportunity to investigate several important phenomena, such as the baryon circulation on 
various scales (from galaxies to clusters), the integral past supernova (SN) activity along 
with the relative role of the two major SN types, the efficiency of gas to galaxies, stars, and 
metals conversion at a cluster scale, etc. (e.g., Vigroux 1977; Matteucci & Vettolani 1988; 
Arnaud et al. 1992; Renzini et ah 1993, hereafter RCDP; Loewenstein & Mushotzky 1996). 
Some of these issues have been extensively discussed in RCDP (see also Renzini 1994) on 
the basis of empirical evidence coming almost exclusively from Einstein X-ray observations. 
In more recent years a great deal of relevant new data have become available from the 
ROSAT and ASCA satellites, and in this paper some of the issues discussed by RCDP are 
revisited, and some new inferences are drawn. 

Clusters are the largest objects on which chemical enrichment can be thoroughly 
studied, and offer the additional advantage of being perhaps the best example in nature 
for which the closed box approximation may hold true. Pei & Fall (1995) have recently 
modeled the chemical evolution of damped Lya absorbers, thus predicting the evolution of 
the global star formation rate with redshift (cosmic time). Global star formation rates all 
the way to z ^ 1 and beyond have been empirically constructed by Madau et al. (1996) 
using data from the Canada-France Redshift Survey (Lilly et al. 1996), and found it in 
remarkable agreement with the predictions of Pei & Fall. Yet, absorbers represent a minor 
fraction of the baryonic matter, unlikely to be well approximated by the closed box model. 
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and Pei & Fall allow for both inflow and outflow in their chemical evolution model. Clusters 
of galaxies are to some extent complementary to Lya absorbers, in that they contain a large 
amount of heavy elements, partly in the hot ICM, partly locked into stars, none of which 
participates in absorbing quasar light. With a bottom-up approach, in this paper the iron 
(and metal) content of clusters and groups is used as a tracer of baryon circulation, past 
star formation, and supernova enrichment at the cluster scale, and beyond. 

The paper is structured as follows: Section 2 presents a recollection of literature data 
concerning the iron and gas content of X-ray clusters and groups, that will form the factual 
basis for the considerations to be developed in the subsequent sections. In Section 3 several 
inferences are derived from the uniformity of these cluster properties, including the evidence 
for clusters having experienced little baryon-exchange with their surrondings, having turned 
gas into galaxies with nearly constant efficiency, and having experienced very similar star 
formation histories. Groups instead present a completely different scenario, which suggests 
that present day clusters cannot have formed by coalescence of groups similar to the present 
day ones. Section 4 deals with the elemental relative abundances at the cluster scale, 
while Section 5 deals with the global metallicity of clusters versus that of the present day 
universe as a whole, and with the global star formation history of the universe. The main 
conclusions of this paper are hnally listed in Section 6. 


2. The Iron Mass to Light Ratio of Clusters and Groups 

Ciotti et ah (1991) introduced the concept of ICM iron mass to light ratio (FeM/L) 
as the ratio of the total iron mass in the ICM over the total optical luminosity 

of galaxies in the cluster, and found values in the range (0.7 — 1.6) x IO^^Mq/Lq for the 
nearby clusters Virgo, Coma, and Perseus. On a wider database, Arnaud et al. (1992) 
found oc Lb, and RCDP hnally adopted Mpg^/Le = 0.01 — O.O 2 M 0 /L 0 as typical 
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for rich clusters. This value has a moderate dependence on the Hubble constant, being 
proportional to (Renzini 1994); for the data presented in this Section h = 1/2 has 
been adopted. The FeM /L relates two pieces of fossil information; the integral amount of 
iron ejected by galaxies in the course of their evolution, and the present luminosity of the 
old stellar remnant of the population that preasumably produced the observed iron, early 
in its evolution. 

The iron mass in the ICM is obtained as the product Zf^y^Micyi of the iron abundance 
times the mass of the ICM gas. The iron abundance is generally estimated at the center 
of the cluster, and one assumes the ICM to be chemically homogeneous. This is the case 
for rich, high temperature clusters, while radial gradients in iron abundance have been 
reported for cooler clusters from ASCA data (Ohashi et al. 1995). In principle, it is 
possible to obtain integrating radially pdV. However, ASCAA PSF does not 

allow a good reconstruction of the cluster ICM density profile. Such integration could be 
accomplished by combining ASCA data with those of other X-ray telescopes with better 
PSF. 

Figure 1 shows Mpg^/Le as a function of the cluster total luminosity. For 
Tb^4 X lO^^L© the FeM/L appears indeed to be constant, with very small scatter around 
an average of ~ 0.02 Mq/L©. However, for Lb~4 x lO^^L© - i.e., in poor groups rather 
than rich clusters - the FeM/L exhibits much smaller values, with no good correlation 
with the cluster optical luminosity. Figure 2 shows the same /L-q data now plotted 

as a function of the temperature of the X-ray ICM. It is now apparent that the FeM/L 
is virtually constant all the way to temperatures down to ~ 1 keV, and then drops 
precipitously by almost three orders of magnitude as the temperature decreases below 1 
keV. The separation between groups and clusters is striking in this diagram. It is worth 
mentioning that groups shown in this figures represent a biased sample in that they are 
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selected for being detected in X rays. Typically, spiral rich groups (similar to the Local 
Group) are not detected in X rays, which may be due to them having an even lower ICM 
content, lower ICM mass to light ratio, and therefore lower FeM/L. 

The drop of the derived FeM/L in poor clusters and groups can be traced to a drop 
in both factors entering in its dehnition, i.e., in the iron abundance and in the ICM mass 
to light ratio. Figure 3 shows indeed ICM mass to light ratio Mjcm/Tb as a function of 
the cluster optical luminosity. While rich clusters have fairly constant ICM mass to light 
ratio, much lower values with large dispersion are exhibited by poor clusters and groups 
with Tb^4 X IO^^Lq. The ICM mass is typically measured within a radius that scales with 
the optical radius of the cluster, and therefore it is smaller for groups than for clusters. 
Figure 4 shows Miqm/Lb as a function of the ICM temperature. Again, fairly constant 
values {Micm/Lb — SOMq/Lq) are exhibited by clusters hotter than ~ 1 keV, while the 
ICM mass to light ratio drops precipitously below ~ 1 keV. 

Figure 5 shows the iron abundance as a function of the cluster optical luminosity. 
Again, rich clusters with Tb^4 x IO^^Lq exhibit a fairly constant iron abundance, at the 
level of ~ 0.3 solar. Poorer clusters instead show a large dispersion, with very low values 
being reached among poorest clusters. The iron abundance is hnally shown as a function 
of ICM temperature in Figure 6. Clusters hotter than ~ 2.5 keV have similar abundances, 
with very small dispersion. Going to cooler clusters Zf^y^ appears hrst to increase up to near 
solar abundance for kT slightly in excess of 1 keV, and then drops precipitously to almost 
zero below 1 keV, with a fairly strong abundance-temperature correlation (see also Arimoto 
et ah 1997, hereafter AMIOR). The presence of this correlation shows that just larger 
statistical errors at low kT cannot account for the apparent large spread of abundances. 
As extensively discussed by AMIOR, very low iron abundances are also derived from X-ray 
data for the ISM of ellipticals, which all have temperatures keV. Also among ellipticals 
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the iron abundance appears to correlate with temperature, the lower the temperature the 
lower the estimated abundance (e.g., Davis & White 1996). 

It is apparent from all these hgures that the FeM/L, the ICM mass to light ratio, 
as well as the iron abundance all correlate much stronger with cluster temperature rather 
than with optical luminosity. Part of the problem may be due to observational errors, as 
cluster optical luminosities are difficult to determine. However, there may be other effects 
at work to make so strong the correlation with temperature. Some will be explored in the 
next sections. 


3. The Baryon Circulation in Clusters and Groups 
3.1. Rich Clusters 

The constancy of the FeM /L and of abundance among clusters give support to the 
notion that galaxy clusters neither lost baryons to, nor acquired baryons from the outside 
in the course of their evolution. Were clusters loosing gas at some stage (e.g., as a result of 
the ICM heating by galactic winds or AGNs) then iron would be lost as well, and sizable 
cluster to cluster differences in the FeM/L should arise. Seemingly, were clusters to accrete 
pristine gas ~ 0) in substantial amount then iron would be diluted, and cluster to 
cluster differences in would arise. Of course, one may argue that all clusters had 

nearly the same amount of baryon exchange with their surrondings, which however seems a 
rather contrived requirement. This latter empirical result is in qualitative agreement with 
the theoretical prediction according to which the baryon fraction of clusters cannot change 
appreciably in the course of their evolution (White et al. 1993; Evrard 1997). 

The constancy of the FeM/L, ICM mass to light ratio, and abundance altogether 
suggest that the conversion of baryons into galaxies, stars and hnally metals took place 



with nearly the same efficiency and at the same epoch in all clusters. Indeed, the constancy 
of the FeM/L says that the same amount of iron was ejected by galaxies per unit present 
stellar population. The most direct interpretation is that there are no appreciable cluster 
to cluster differences in the global stellar initial mass function (IMF) and star formation 
history, unless variations in IMF are rather precisely compensated by variations in the 
distribution of stellar ages (again, a rather contrived requirement). For example, a difference 
by Ax = 0.5 in the IMF slope would imply a factor ~ 4 difference in the FeM/L (cf. 
RCDP). Differences from cluster to cluster in the average age of the stars producing the 
present optical luminosity of the cluster would also affect the FeM/L. Assuming the cluster 
luminosity evolution to be dominated by the passive aging of the elliptical galaxy population 
{L oc~ factor of two difference in average stellar age implies a similar difference 

in luminosity, hence in the FeM/L. 

As well known, most of the cluster light Lb is provided by elliptical galaxies and 
galactic spheroids. Moreover, they represent an even larger fraction of the stellar mass 
content of clusters since these passively evolving systems have larger M/ Lb ratios than star 
forming galaxies. It is now well established by several independent lines of evidence that 
the bulk of stars in cluster ellipticals formed at high redshift, i.e., 2:^3 (for a review see 
e.g., Renzini 1995). This comes from the tightness of the color—a relation of ellipticals in 
nearby clusters (Bower, Lucey, & Ellis 1992), the tightness of the distributions of elliptical 
galaxies about their fundamental plane at low and high redshift (Renzini & Ciotti 1993; 
van Dokkum & Franx 1996; Pahre, Djorgovski, & de Carvalho 1997), the tightness of the 
Mg 2 — a relation, again at low as well as high redshift (Bender, Ziegler, & Bruzual 1996), 
the existence of red (elliptical) galaxies at high redshift (Hamilton 1985; Aragon-Salamanca 
et al. 1993; Dickinson 1995), and the tightness of the color-magnitude relation of cluster 
ellipticals at high redshift (Ellis et al. 1996). The possibility has been advocated of hue 
tuning between age and metallicity effects conspiring to preserve the tightness of these 
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relations at low-z while allowing for a large age spread (Worthey, Trager, & Faber 1995). 
This appears to be rnled ont by the observation that the tightness is preserved even for 
clusters at large lookback times, while it would obviously be destroyed if a large age spread 
was present among ellipticals in low-;^ clusters (Kodama & Arimoto, 1997). 

It is most straigtforward to identify in the spheroidal populations the producers of 
the bulk of the ICM iron in clusters, since most of rich cluster stellar luminosity comes 
from ellipticals, SOs, and the bulges of early-type spirals (which are also dominated by very 
old stellar populations, cf. Ortolani et ah 1995; Jablonka & Alloin 1995). If so, then the 
bulk of iron had to be produced by stars and ejected from galaxies at very early times, 
say at z^2 — 3, even if a major fraction of iron could have been produced by late SNIa 
exploders (RCDP). This - coupled with a non-evolving baryon fraction of clusters (White 
et al. 1993) ~ implies an ICM iron abundance independent of redshift, all the way to very 
high redshifts. Some evidence in this direction is indeed emerging, though for the moment 
limited to moderate redshifts (cf. Figure 6). 

Although it appears most natural to attribute to ellipticals and spheroids the prime role 
as iron (metal) producers in clusters, yet other alternatives have been entertained and are 
worth considering. Giant elliptical galaxies and large spheroids have a rather deep potential 
well, hence ejecting large amounts of metals requires a great deal of energy. Dwarfs instead 
de-gas much more easily than giants, having much shallower potential wells, and therefore 
dwarfs (e.g., dSph galaxies) are likely to have ejected more iron per unit present light 
than giants. This is indeed a common feature of galactic wind models for the formation 
of ellipticals and spheroids (e.g., Larson 1974; Arimoto & Yoshii 1987; Matteucci 1992). 
Clearly, the problem is whether there were enough dwarfs to make an appreciable fraction 
of the ICM iron. It has actually been argued that most ICM itself might have been ejected 
from dwarfs, given that the dwarfest among dwarf spheroidals have a very low baryonic 
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fraction, hence had to eject most of their baryons if they were to start with the cosmic 
share (Trentham 1994). However, both the faint end of the galaxy luminosity function, and 
the relative amounts of baryons and metals ejected by dwarfs are sufficiently uncertain to 
make very hard either to prove or disprove this scenario. Most baryons are now in the ICM 
rather than within galaxies (i.e., galaxies contain ~ 20% of the baryons, e.g., David et ah 
1990; White et ah 1993). Hence, the question is whether the baryons were separated from 
their (small scale) dark matter complement by local starbursts, or by mutual stripping in 
the process of hierarchical clustering (Nath & Chiba 1995), when small scale density peaks 
dissolved inside the collapsing cluster. In the former case a population of low metallicity 
stars would be present, either still clustered in e.g., dSph’s, or dispersed through the whole 
cluster potential well. However, a dispersed population can also arise from tidal stripping of 
galaxies in the course of the cluster evolution (cluster harrassment, cf. Moore et al. 1996), 
and it will be very hard to distinguish observationally between the two options. Upper 
limits to a diffused stellar population in clusters have been reported to be at the level of 
~ 20% of the total cluster light (Melnick, White, & Hoessel 1977). A promising way to 
better determine this fraction may come from detecting and counting isolated planetary 
nebulae that are now being found in clusters (Theuns & Warren 1997; Arnaboldi et al. 
1996). Theuns & Warren argue that perhaps as much as 40% of the cluster stellar light 
comes from a dispersed population, but also point out that this estimate is based on an 
uncertain frequency of planetary nebulae per unit light of the parent population. This is 
known to vary by a very large amount among spheroidals, with most metal rich ones having 
a several times smaller PN productivity compared to metal poor ones (Hui et al. 1993; 
Ferguson & Davidsen 1993). 

Another possibility is offered by an early generation of massive stars in a now extinct 
(i.e.. Lb — 0) population produced in a bimodal star formation scenario (Arnaud et al. 
1992; Elbaz, Arnaud, & Vangioni-Flam 1995). One postulates that only massive stars (e.g.. 
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M>2 Mq) would form in starbursts, while the low mass stars that shine today in galaxy 
clusters would come only later in a quiescent star formation mode. For this model to work 
hne tuning is required to explain the constancy of the FeM/L, with the ratio of the star 
yields of the two star formation modes being the same in all clusters. Moreover, there 
is little evidence for bimodal star formation in nature. For example, a typical globular 
cluster formed ~ IO^Mq of stars in probably less than 10 ® yr, with a star formation rate 
of ~ lM 0 yr“^ per cubic parsec (!), a strong starburst indeed, and yet they were able 
to produce the wealth of low mass stars without which we would not see them today. 
Moreover, abundant low mass stars have been detected in Galactic regions with active 
high mass star formation (Zinnecker, McCaughrean, & Wilking 1993), and therefore the 
occurrence of bimodal star formation appears rather implausible. All in all, it appears 
natural to adopt the view that the metals we see in clusters were produced by the same 
stellar population which low mass component allows us to see cluster galaxies today. 

3.2. Galcixy Groups 

An inspection to Figures 1-6 immediately reveals that what holds for clusters, either 
hotter than ~ 2 keV or brighter than ~ 4 x IO^^Lq, apparently does not hold for 
cooler/fainter clusters and groups. I hrst assume and Mjcm/Fb ratios for these latter 
objects at face value, and follow their astrophysical implications. Some of such implications 
were already mentioned in RCDP, even if at that time data were available for only one 
object, namely the NGC 2300 group (Mulchaey et ah 1992). The very low values of 
the FeM/L shown in Figure 2 imply that a great deal of iron should have been lost by 
the groups, hence a major fraction of their baryons along with it. This appears to be in 
agreement with the IGM mass to light ratio being indeed much lower than in clusters 
(Figure 4). However, this cannot be the whole story, because besides a lower FeM/L several 
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groups appear to have a much lower iron abundances compared to clusters (Figure 6). It 
seems as if, after having shed a major fraction of their original baryon content, such groups 
had to re-accrete a modest amount of baryons from their surrondings, where pristine gas 
largely diluted the metal contaminated ejecta of the groups themselves (RCDP). This 
alternate flow of baryons, out and in groups, may arise as a result of a declining galactic 
wind activity inside groups, and at least at early times such activity could well have been 
strong enough to drive gas out of groups (Renzini 1994). However, the whole scenario 
appears rather contrived, as it may be difficult for groups to reaccrete gas once its adiabat 
has been raised by energy input from the galactic winds (Kaiser 1991). 

There is however another interpretation still open. It is apparent from Figures 2, 4 
and 6 that for kT>2 keV all plotted cluster properties are nearly constant and exhibit very 
small scatter. Cooler clusters and groups, especially for kT<l keV, show instead a wild 
range of abundances, gas mass to light ratios, and therefore FeM/L values. It is worth 
emphasizing that the diagnostic tools to get the iron abundance are radically different in the 
hot clusters compared to the cooler ones. Above kT ~ 2 keV iron abundances are derived 
primarily from the iron-K complex at ~ 7 keV, where lines arise from transitions down to 
the K shell of He-like and H-like iron ions. The atomic physics used in this plasma emission 
models is therefore rather simple. At cooler temperatures instead, iron abundances are 
derived from the iron-L complex at ~ 1 keV, where lines arise from transitions down to the 
L shell of progressively more and more complex iron ions as temperature becomes lower and 
lower. The more complex the ions, the less secnre the involved atomic physics calculations, 
especially the collisional excitation probabilities, and it appears legitimate to entertain the 
suspicion that much of the observed trends below ~ 2 keV may just be an artifact of some 
systematic error in the iron-L diagnostics (cf. AMIOR for an extensive discussion). In the 
case of the Virgo cluster [kT ~ 2.9 keV) it has been shown that the iron abundance derived 
from iron-L lines is in very good agreement with that derived from iron-K lines (AMIOR). 
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Hwang et al. (1997) come to the same conclusion analysing a sample of clusters with 
kT in the 2-4 keV temperature range. Unfortunately, by no means the reliability of the 
iron-L diagnostics at kT ~ 3 keV can prove its reliability at lower temperatures, since then 
other, much more complex iron ions are involved (AMIOR). In conclusion, the drop of iron 
abundance below ~ 1 keV may not be real at all, and therefore it is premature to conclude 
that groups had to re-accrete baryons at late epochs. It appears that the best way to assess 
the reliability of iron-L diagnostics at low temperatures will be offered by collecting ASCA 
observations of objects for which the iron abundance is independently known, such as stars, 
supernova remnants, starburst galaxies etc. (cf. AMIOR for a preliminary discussion). 

While the iron abundances may be in error at low kT values, it is worth emphasizing 
that the low Micm/Lb ratios shown in Figure 4 should be only marginally affected by 
possible systematic errors in the iron-L diagnostics, i.e., while abundances and FeM/L 
values may be severely underestimated below ~ 1 keV, the low values of the ICM mass to 
light ratios in groups are more robust, and one can safely conclude that groups are genuinely 
gas poor, likely to have lost a major part of their original share of baryons. Alternatively, 
in groups baryon conversion to stars would have been much more efficient than in clusters, 
proceeding almost to the last “drop” of them, which seems a less plausible interpretation. 


3.3. Clusters vs. Groups and the Formation Epoch of Clusters 

The systematic difference in the Miqm/Lb ratio of groups relative to rich clusters can 
have far reaching consequences. Indeed, for merging groups the Micm/Tr ratio should stay 
the same or decrease slightly (Evrard 1997), therefore failing to meet the high Micm/Tb 
values typical of clusters, which are up to ~ 30 times higher than those in groups. This 
conclusion would be further reinforced were the iron-L diagnostics correct, since merged 
groups would also have much lower FeM /L compared to clusters. 
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The inference is that rich clusters cannot have formed by merging groups similar to 
those in the nearby universe. If clusters formed by hierarchical merging, merging should 
have taken place when groups were still gas rich, likely at high redshift. This favors low-f2 
universes, where little merging and cluster growth takes place at late epochs (low z), while 
most of merging activity “switches off” at 1 + ~ ~ 5 in an open universe, or at 

1 + 2 : ~ ~ 1.5 in a ACDM universe (White 1997). Late accretion of groups onto rich 

clusters may still take place, but should occur at such a low rate to leave substantially 
unaffected the ICM mass to light ratio of the clusters. Independent evidence for very 
little or no cluster evolution all the way out to 2 ~ 0.6 is now emerging both from optical 
(Carlberg et al. 1996) as well as X-ray selected cluster samples (Rosati et ah 1997). 

This conclusion is further reinforced by some of the considerations in Section 3.1. 
Indeed, with most of stars in clusters having formed at high redshift ( 2 ^ 8 ), and with the 
concomitant production and ejection of iron, much of the heating of the ICM also took 
place at such early epoch (e.g., Renzini 1994). If clusters were not already in place such 
heating would have drasticaly affected the baryon content of clusters, as it appears indeed 
to have done for groups. While small cluster to cluster variations in the baryonic fraction 
might have been detected (Loewenstein & Mushotzky 1996b), such variations are tiny 
compared to the range covered by the Miqm/Lb ratio of groups in Figures 3 and 4. Rich 
clusters appear to be “older” than the bulk of the old stellar population by which they are 
dominated. 

However, one may also envisage a scenario in which groups destined to assemble 
forming a rich cluster expelled gas at an early stage, the gas remained conhned as an 
intragroup medium, and hnally was collected inside the cluster along with the groups. 
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4. The Elemental Ratios in the Intracluster Medium 

It is now well established that in the metal poor stars of the Galactic halo the 
a-elements are enhanced with respect to iron relative to the solar proportions (e.g., 
Wheeler, Sneden, & Truran 1989; Bessell, Sutherland, & Ruan 1991), and a similar 
enhancement is observed in the metal rich stars of the Galactic bulge (McWilliam & Rich 
1994). The current interpretation of this o-element overabundance in the whole Galactic 
spheroid appeals to the prompt release of a-elements by short living massive stars producing 
Type 11 SNs, coupled to the somewhat delayed release of a major fraction of the whole iron 
by Type la SNs (e.g., Greggio & Renzini 1983; Matteucci & Greggio 1986; Ruiz-Lapuente, 
Burkert, & Ganal 1996). Various degrees of a-element enhancement can thus be obtained, 
depending on the adopted time scale of the iron release by SNla’s relative to the time scale 
of star formation. This scenario is referred to as the “standard chemical model” for the 
Galactic chemical evolution (RGDP). Moreover, population synthesis methods indicate that 
an a-element enhancement may also be present in ellipticals, reaching up to 0.2-0.3 dex in 
the most massive ones (e.g., Worthey, Faber, & Gonzales 1992; Davies, Sadler, & Peletier 
1993). By analogy with the Galactic spheroid, the favored interpretation has been in terms 
of the relative roles of the two SN types during the fast formation process of ellipticals. 

A fast completion of star formation (on time scales of, say ^ few 10® yr) implies that 
a fraction of the iron released by SNla’s should flow directly out of galaxies during an 
early galactic wind phase, without ever being incorporated into stars. If so, a chemical 
asymmetry should be established between galaxies and the IGM, with galaxies being slightly 
overabundant in a-elements relative to iron, and the IGM being slightly overabundant 
in iron relative to the a-elements (RGDP). The size of this asymmetry is hard to predict 
theoretically, as it depends on the relative time scales of star formation on the one hand, 
and of the SNla iron release on the other, with both time scales being poorly known. The 
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effect is not expected to be large. If all the iron from SNIa were ejected from galaxies the 
a-dehciency in the ICM would be [Q;/Fe]~ —0.4 (RCDP). The actual dehciency may be 
much smaller than this, since part of the SNIa iron is likely to be incorporated into stars, 
while SNII products are also ejected from galaxies. 

Early attempts at measuring the abundance of a-elements in the ICM indicated a 
possible overabundance of oxygen relative to iron (Canizares et ah 1982) and a near 
solar Si/Fe ratio (Mushotzky et ah 1981), but both with large statistical and systematic 
uncertainty. ASCA has now provided much better data, and Mushotzky (1994) has initially 
reported a fairly high a-element enhancement, with < [Q;/Fe]>~ +0.4 (from his Table 
3). More recently Mushotzky et al. (1996) have revised down this estimate, and report 
detailed elemental abundances for several kT =3-4 keV clusters, all showing a moderate 
a-element enhancement. Uncertainties in the abundances are reported to be of about a 
factor of two, at the 90% conhdence level. Taking a global average for O, Ne, Mg, Si, 
and Fe, one obtains a modest < [Q;/Fe]>~ +0.2. This is still in conflit with the predicted 
asymmetry, and argues for not only the galaxies, but the ICM as well being dominated 
by SNII products (Loewenstein & Mushotzky 1996a). If so, this would demonstrate the 
impossibility to extend to galaxy clusters the standard chemical model that apparently 
holds for the Galactic nucleosynthesis (AMIOR). 

However, Ishimaru & Arimoto (1997) have recently pointed out that the small 
a-element enhancement in the ICM comes from Mushotzky et al. (1996) having assumed 
reference solar abundances from “photospheric” model atmosphere analysis. The result is 
different if one uses “meteoritic” abundances instead, because the meteoritic iron abundance 
is ~ 0.16 dex lower than the photospheric value (Anders & Grevesse 1989). After noting 
that the meteoritic iron is now more generally adopted, Ishimaru & Arimoto conclude that 
the a-element enhancement in the ICM virtually disappears, going down from ~0.2 to only 
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~ 0.04 (consistent with no enhancement or even small depletion, given the reported errors). 

At this stage it is safe to conclnde that there is no strong evidence for an enhanced 
a-element proportion in the ICM, and that the galaxy-ICM chemical asymmetry - if it 
exists - must be small and probably hard to detect given the current errors in abundance 
determinations. Abundance ratios in the cluster as a whole (including stars and ICM 
together) are therefore very close to solar, as expected in the frame of the standard chemical 
model (RCDP). Therefore, the applicability of the standard chemical model to clusters 
is not invalidated, and the possibility of a major contribution from SNIa’s to the iron 
enrichment of clusters remains viable. 


5. The Metallicity of the Present Day Universe and the Past History of Star 

Formation 

Clusters of galaxies are the largest entities for which we can measure the metallicity. 

It is important to notice that the mass of iron in the ICM is comparable to that locked 
into galaxies, having assumed the average abundance of stars in galaxies to be solar 
(RCDP). While the abundances in either the ICM or galaxies are independent of the 
distance scale, the cluster average abundance does actually depend on the assumed Hubble 
constant, because so do both the mass in galaxies and Micm, and each of them in a 
different way. I assume as prototypical the Coma cluster values adopted by White et ah 
(1993): Micm — 5.5 x and M* ~ and derive for the cluster iron 

abundance: 

..P. ZSmVicm + ZfM, b.bZ^i^h-^1^ + 

“ A/ICM + A^. 5.5A-5/2 + ft-i ’ ' ’ 

where Zf® is the average abundance of stars in galaxies and M* is the mass in stars. 

With Zf^-^ = 0.3 solar and Zf® = 1 solar, equation (1) gives a global cluster abundance 

of 0.34, 0.37, and 0.41 times solar, respectively for h = 0.5, 0.75, and 1. Under the same 
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assumptions, the ratio of the iron mass in the ICM to the iron mass locked into stars is: 

Z^‘M. 

or 4.6, 2.5, and 1.65, respectively for h = 0.5, 0.75, and 1. Note that with the adopted values 
for the quantities in equation (2) most of the iron is in the ICM, rather than now locked 
into stars, especially for low values of Ho. These estimates could be somewhat decreased if 
clusters contain a sizable population of stars not bound to censed individual galaxies (cf. 
Section 3.1). The iron content of galaxies may also have been underestimated, because 
so does a luminosity-weighted abundance compared to the mass-weighted abundance 
(e.g., Greggio 1997). In addition, the galaxy (baryonic) M^/L ratio might be higher than 
estimated by White et al. (1993), i.e., <M^/L-q >= 6.4/i. All these effects together may 
contribute to reduce the ICM to galaxies iron ratio below the somewhat embarrasingly 
large values given above. 


With the adopted masses and iron abundances for the two baryonic components one 
can also evaluate the total cluster FeM/L. One needs to specify the average M*/Lb ratio, 
and adopting for consistency the value given by White et al. (1993), one gets: 

i\/rlCM I A/T* 

^ ~ 1.3 X 10-2(1.65 (Mq/Lq), (3) 

Lb 

or FeM/L=0.037 or 0.034 Mq/Lq, respectively for h = 0.5 and 1. The total FeM/L is 
therefore fairly insensitive to the adopted distance scale, and is close to the value previously 
estimated (i.e., O.OOMq/Lq, RCDP; Renzini 1994). Simple calculations (cf. RCDP) show 
that to reproduce this value one needs either a fairly flat IMF [x ~ 0.9) if all iron is 
attributed to SNII’s, or a major contribution from SNIa’s, if one adopts a Salpeter IMF 
{x = 1.35). The former option dictates a substantial a-element enhancement, similar to 
the values observed in the Galactic halo ([Q;/Fe]~ -1-0.5). The latter option instead predicts 
near solar proportions for the cluster as a whole. On the basis of the discussion in Section 
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4 one concludes in favor of the second option. Moreover, the total metal mass to light ratio 
of clusters is ~ 0.3 Mq/Lq, given that in the solar proportion iron accounts for about 10% 
of all metals. There is no doubt that the bulk of such metals other than iron - hence of the 
metals as whole - was produced by SNIIs. Now, the number of SNIIs exploded at early 
times in a stellar population that has faded to luminosity Lb when aged to ~ 15 Gyr is 
~ 0.1 X Lb/Lq (RCDP). Therefore, to produce the observed metal mass to light ratio of 
clusters each SNII must have contributed on average ~ 3Mq of heavy elements, which is in 
agreement with current nucleosynthesis calculations (e.g., Woosley & Weaver 1995). 

A critical issue is to what extent the cluster global metallicity, and the ICM to galaxies 
iron share are representative of the low—universe as a whole. For example, Madau et 
ah (1996) adopt Ho = 50, a stellar mass density parameter G* = 0.0036, a baryon mass 
density parameter Gb = 0.05, an average solar metallicity for the stars, and a negligible 
metal content for the intergalactic medium (IGM), that comprises the vast majority 
of the baryons. With these assumptions the metallicity of the present day universe is 
~ 1 X 0.0036/0.05 = 0.07 solar, or ~ 5 times lower than the measured value in clusters of 
galaxies. In the same frame, the fraction of baryons in galaxies (stars) is also ~ 7% globally, 
which compares to 1/(1 + 5.5/1-3/2) in clusters, or ~ 6% and ~ 10%, respectively for 
h = 0.5 and 0.75. Therefore, it appears that the efficiency of baryon conversion into galaxies 
and stars (G^/Gb) adopted by Madau et ah (1996) is nearly the same as that observed in 
clusters, which supports the notion of clusters being representative of the low-;^ universe. 
The metallicity of the clusters is however ~ 5 times higher than the metallicity of the low-;^ 
universe as conservatively adopted by Madau et ah. The difference comes from having 
attributed to the IGM a very low metallicity, hence assuming held galaxies losing only a 
negligible amount of metals. This implies a factor ~ 5 lower efficiency in metal production 
per unit mass turned into stars. I will discuss the two aspects in turn. 
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Is there any reason why most of the produced metals should be ejected by cluster 
galaxies, and instead fully retained by held galaxies, while reaching the same average stellar 
metallicity? The only hint for a systematic cluster/held difference may be offered by ram 
pressure stripping, that to some extent should be active in clusters but not in the held. 
However, in clusters iron was most likely ejected from galaxies as a result of supernova 
heating, rather than stripped by ram pressure (cf. RCDP). In particular, ram pressure 
stripping cannot have played a major role because as shown by Figure 2 the FeM/L is 
the same in clusters with moderate velocity dispersion (or equivalently, kT ~ 2 keV) as in 
clusters with high velocity dispersion (or kT ~ 10 keV). Were ram pressure important one 
would have expected the FeM/L to increase with ICM temperature. There appears also to 
be no strong argument for held galaxies having shed much less iron than cluster galaxies. 
The small starburst galaxies making the excess blue galaxy counts (Lilly et al. 1995) are 
indeed likely to have ejected a sizable amount of metals before “fading to oblivion”. 

Seemingly, there are no strong arguments supporting the factor of ~ 5 diherence 
in metal production efficiency, unless one is willing to postulate a hatter stellar IMF in 
cluster relative to held galaxies. The diherence might be somewhat smaller if one allows 
for the possibility that there are more stars in clusters than adopted here (cf. Section 
3.1). However, if such stars escaped detection within well studied clusters, they may have 
gone unnoticed in the held as well. Therefore, it would appear quite ad hoc to appeal to a 
much higher metal productivity of stellar populations in galaxy clusters compared to their 
held counterpart, and it appears reasonable to conclude that the global metallicity of the 
present day universe may well be nearly the same as that observed in galaxy clusters, i.e., 
~ 0.3 — 0.4 times solar. If so, there should be a comparable share of metals in the held 
IGM, as there is in the cluster ICM. 


While it appears that the global, time-averaged rate of metal production in the universe 
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may have been underestimated by perhaps as much as a factor ~ 5, this does not necessarily 
imply that the global, time-averaged rate of star formation was also underestimated by 
the same factor. More simply, the metal productivity of galaxies per unit of their present 
mass may have been underestimated by Madau et ah, if the cluster galaxies productivity 
is to be taken as representative of the universe as a whole. One cannot exclude that actual 
average rate of star formation may have been underestimated, which would be the case if a 
population exists of unaccounted stars in clusters and in the held. 


6. Conclusions 

In this paper several relations coming from combining galaxy cluster X-Ray and optical 
data have been presented and discussed, and astrophysical inferences have been drawn from 
them. The main results can be summarized as follows. 

1. The ICM iron mass to optical light ratio of clusters with kT>2 keV appears to be 
constant, with a small dispersion which is fully consistent with observational errors. Among 
these clusters also the iron abundance and the ICM mass to light ratio appear to be 
constant with small dispersion. 

2. Among poor clusters and groups, instead, the FeM/L, the ICM mass to light ratio, and 
the iron abundance all drop by orders of magnitude compared to the corresponding values 
typical of the rich clusters. This is especially evident for ICM temperatures below ~ 1 keV. 

3. The constancy of the FeM/L among clusters indicates that clusters did not lose an 
appreciable fraction of their original share of the cosmic baryons. It also argues for the 
stellar initial mass function being nearly the same in all clusters, as well as for the average 
stellar ages being also the same from cluster to cluster, within less than a factor ~ 2. 


4. As a major fraction of the cluster optical luminosity comes from ellipticals and bulges. 
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the stellar population responsible for the production of the iron and the other heavy 
elements now in the ICM is most naturally identihed with the high mass tail of the same 
population now surviving in ellipticals and bulges. Since such spheroid populations formed 
at high redshift much of the iron should have been ejected from galaxies at a very 

early stage. This is in agreement with the observed iron abundance in moderate redshift 
clusters being the same as in local clusters. 

5. The possibility remains that a non-negligible fraction of the cluster stellar population 
resides out of censed galaxies. If so, a corresponding fraction of the ICM metals could 
have been manufacted by this dispersed stellar population. In principle, part of the metals 
could have been produced by a now extinct stellar population, generated in a bimodal star 
formation mode. It is argued that this hypothesis is rather implausible. 

6. The drop of the ICM mass to light ratio in groups indicates that - contrary to clusters - 
they have experienced a major loss of gas (baryons) in the course of their evolution. The 
very low metallicity indicated by current X-ray diagnostics for most of these groups, if real, 
suggests that baryons should have been re-accreted after an initial de-gassing, which seems 
rather contrived. However, it is possible that current iron-L diagnostics for kT<l keV is 
affected by systematic errors. 

7. The large differences between the cluster and the group ICM mass to light ratios suggest 
that clusters did not form by agglomerating groups similar to the present day ones. If such 
agglomeration took place, it must have occurred while the groups were still gas rich, i.e., 
before major star formation took place thus propelling part of the baryons out of them. 
This argues for the assembly of rich clusters having been completed at high redshift, which 
favors low-H universes. 

8. There is not much evidence for a predicted chemical asymmetry in the a-element to iron 
ratios between the ICM and the stellar component of clusters. However, data are consistent 
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with no a-element enhancement in the ICM, and therefore with a combined SNIa and SNII 
enrichment to give near solar elemental proportions on the whole cluster scale, as predicted 
by the standard model for the Galactic chemical evolution. In particular, there appears to 
be no need to invoke a special IMF or a suppression of SNIa’s in clusters, compared to our 
own Galaxy. 

9. It is emphasized that the overall cluster metallicity is about 5 times higher than the 
currently adopted average metallicity of the present-day universe, in spite of a similar 
fraction of the baryons having been converted into galaxies and stars in the clusters as well 
as in the general held. Possible origins of the difference (or discrepancy) are discussed, 
including an underestimate of the efficiency of metal production per unit mass of baryons 
converted into stars. 

10. The amount of iron (metals) in an undetected, probably hot intergalactic medium 
comprising most of the bayons is predicted to be comparable to - or even larger than - the 
amount of iron now locked into stars in the present day universe. The metal abundance of 
this IGM is predicted to be ~ 1/3 solar, if the metal productivity of stellar populations is 
the same in the held as it is in clusters. 

I am grateful to Mauro Giavalisco and Piero Madau for their support and for useful 
discussions about the metallicity of the present day universe and the integrated past 
star formation activity, and to Ralph Bender and an anonymous referee for constructive 
comments. I wish to thank Richard Ellis and Trevor Ponman for stimulating discussions, 
and for their invitation to participate at the Royal Astronomical Society meeting of October 

11. 1996, and the Pontihcial Academy of Sciences for its invitation to attend the Vatican 
Gonference on “The Emergence of Structure in the Universe at the Level of Galaxies”, 
November 25-29, 1996, as at these meetings I had the opportunity to report and discuss in 
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Fig. 1.— The iron mass to light ratio of the ICM of clusters and groups as a function of 
the total optical luminosity Lb of the cluster galaxies. Data are taken from the following 
sources: hlled circles: Arnaud et ah (1992); hlled triangles: Tsuru (1993); open triangle: 
David et al. (1994a); open square: Mulchaey et al. (1993); hlled square: Ponman et al. 
(1994); open circles: Mulchaey et al. (1996). 

Fig. 2.— The same as Figure 1 but as a function of the ICM temperature. 

Fig. 3.— The ICM mass to light ratio, i.e., the mass of the ICM per unit light of the cluster 
galaxies, as a function of the total optical luminosity Lb of the cluster galaxies. The same 
objects as in Figure 1 are displayed. 

Fig. 4.— The same as Figure 3 but as a function of the ICM temperature. 

Fig. 5.— The iron abundance in the ICM for the objects in Figure 1 as a function of the 
total optical luminosity Lb of the cluster galaxies. 

Fig. 6.— The same as Figure 5 but as a function of the ICM temperature. Data for six 
clusters at moderately high redshift (< 2 : >~ 0.33) are also included and represented by 
small hlled circles (Allen et al. 1996; Donahue 1996; Matsuura et al. 1996; Schindler et al. 
1997). Note that for kT >~ 3 keV the iron abundance is derived from the iron-K complex, 
while for lower temperatures the iron-L complex is used. 
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